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Abstract 

The higher order super symmetric partners of the Schrodinger's Hamil- 
tonians can be exphcitly constructed by iterating a simple finite difference 
equation corresponding to the Backlund transformation. The method can 
completely replace the Crum determinants. Its limiting, differential case of- 
fers some new operational advantages. 
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1 Introduction 



It is notable that while the exact solutions in general relativity fill already ample 
handbooks, the exactly solvable problems of quantum theory form still an excep- 
tional area. Yet, it seems that they are just the "top of an iceberg", an opinion 
recently confirmed by advances in the exact methods [|I|]. One of the simplest ways 
of obtaining new exactly solvable spectral problems in Schrodinger's quantum me- 
chanics is to look for pairs of quantum Hamiltonians coupled "supersymmetrically" 
by an intertwining operator Ai 

H,A, = A,Ho. (1) 

This equation should hold on a certain dense domain of the Hilbert space of states 
7i, assuring that the spectrum of Hq yields some information about the spectrum 
of Hi and vice versa. The best known cases of occur in = L^(R), where Hq 
and Hi are the Schrodinger Hamiltonians 

and Ai, A| are the first order differential operators 

Ai = ^ + A(x), 4 = -^+A(a:). (3) 

By introducing (H) and @ into ([|), and after an elementary integration, one gets 
the Riccati equation for 

-(3[(x)+/3',ix) = 2[Voix)-e], (4) 

where e is an arbitrary integration constant called the factorization energy (a ter- 
minology which does not imply that e must be a physical energy). Simultaneously, 

Vi{x) = Vo{x)+(3'i{x). (5) 

Equations (H) and are necessary and sufficient conditions for the Hamiltonians 
to be factorized as 

HQ-e = ^AlAi, Hi-e = ^AiA\. (6) 

The susy partners (l^)-® appear in all quantum problems solvable by factorization 
[0. An element of arbitrariness in this construction was explored in 1984 [Q; its 
presence in Witten's susy quantum mechanics was noticed by Nieto links with 
the historical Darboux theory |^ were established by Andrianov et al 0. 

The higher order intertwining was considered by Crum ^ (see also Krein 
Adler and Moser p). In physical applications though, for a long time, the atten- 



tion was focused mostly on the first order intertwining operator Ai (see |]T0| and 



2 



references quoted therein). The idea of the higher order Ai reemerged in 80'tieth 



nj |l^] and it was systematically pursued since 1993 (see the works of Veselov and 



Shabat [jT3|, Andrianov et al [|l4l, Eleonsky et al ||T5|, Bagrov and Samsonov [Q, 
Fernandez et al [[iTi [18|, and one of us [|l9l). As it turns out, the higher order 
susy partners of an arbitrary Hq can be generated just by iterating the Darboux 
transformation with different factorization constants e (see interesting theorems by 
Bagrov and Samsonov in [|l6l). Formally, this can be done by introducing a higher 
order intertwiner constructed by means of the determinant formula of Crum |^, 
a method fundamentally important, though not very easy to apply. The purpose 
of this note is to explore a faster method based on the Backlund transformation 
[pO| which operates very simply at the level of the /3-functions (superpotentials) in 
equations (^-(U). It consist just in applying the finite difference calculus, without 
caring about wave functions, differential equations, etc. Moreover, if applied to a 
sufficiently rich initial family of solutions to (|^), it permits to recover the complete 
parametric dependence in all subsequent susy steps without the need of integrating 
any differential equations |]21], ^ . We also show that the differential case of the 
algorithm gives an entire susy hierarchy for a fixed e just in terms of quadratures 
and algebraic operations. 



2 The finite diff^erence algorithm 

The solutions of the Riccati equation (H), and the potentials (^, depend implicitly 
on the factorization energy e G M. We shall make explicit this dependence by 
writing /3i(x,e) for the solutions to (^), and Vi{x,e) for the corresponding susy 
partners of Vo(x). Hence, the intertwiner in (|1|) reads Ai{e) = d/dx + /3i(x,e), 
while the corresponding Hamiltonians are Hq = (l/2)A|(e)Ai(e) + e, and -ffi(e) = 
{l/2)Ai{e)A\{e) + e. Notice that Hq does not depend on e, yet it admits many 
different factorizations for many factorization energies e, each one defining a different 
intertwined Hamiltonian Hi{e), with a different potential Vi{x,e). 

Let us fix now e = ei, and consider one of these potentials: 

V,{x,e,) = VQ{x)+(3[{x,e,). 

Since Vi(a;,ei) is a susy partner of Vq{x), the eigenf unctions ipi of the Hamiltonian 
Hi{ei) can be generated from the eigenf unctions tpQ of Hq by using the operator 
Ai{ei). It means simply that if iljQ{x,e) is a sufficiently smooth function and if 

ipi{x,e) oc Ai{ei)ipQ{x,e), then: 

Hqi/jq^x, e) = eipo{x, e) =^ Hi{ei)ilJi{x, e) = eV'i(x, e). (7) 

Notice that ([7p does not require the square integrability of either ipQ or ipi. Never- 
theless, if e corresponds to one of the discrete spectrum eigenvalues E of Hq, then 
ipQ can be chosen square integrable. 
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A question arises now, can Vi(x, ei) be a convenient point of departure for the 
next susy step? The answer is positive and admits an exphcit construction. In 
geometric form it was described by Adler as a mapping of complex polygons corre- 
sponding to the Backlund transformations |^ (see also Veselov and Shabat |]13|). 
Its presence in the supersymmetric algorithm can be shown directly, using only the 
Riccati equation (H). Indeed, if /5i(x,e) is a solution of (^), the function 



ui[x, e) = exp 



Pi{x, e)dx 



satisfies HoUi{x,e) = eui{x,e). Therefore, the function defined as 



— / /3i{x, e)dx 



U2{x, ei, e) = Ai(ei)Mi(x, e) = ei) - pi{x, e)] exp 



satisfies Hi{ei)u2{x,ei,e) = eu2{x,ei,e), and the function P2{x,ei,e) defined by 



U2{x, ei, e) = exp 



l32{x, ei, e)dx 



must fulfill the new Riccati equation: 

- P'2ix, ei, e) + /3i{x, ei, e) = 2[Vi{x, d) - e]. 
Reading back and using (H) one has 

P2{x,ei,e) = -— lnM2(x,ei,e) = A(x,e) - — . _ . t: 

and using again one ends up with an algebraic expression 

2(ei - 6) 



(9) 



(10) 



/32(a;,ei,e) = -/5i(a;, ei) 



/5i(x,ei) -/3i(x,e)' 



'111 



which is the basic element of the auto-Backlund transformation for the Korteweg- 
deVries (KdV) equation (as reported in see also and p6[). Thus, the /?- 
functions of the second susy step are algebraically determined by a finite difference 
operation performed on the /9-functions of the previous step. They give the next 
operator A2{ei,e) = d/dx + P2{x,ei,e) intertwining the Hamiltonian -ffi(ei) with a 
new one H2{ei, e) = -(1/2) d^ /dx^ + ^2(0;, ei, e). 



H2{ei.e) A2(ei,e) = A2(ei,e)ifi(ei) 



where 



V2{x, ei, e) = Vi{x, ei) + [32{x, ei, e). 
Quite obviously (|1]) and (0) imply 

H2ieue)A2 = A2H0, 



(12) 
(13) 
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where A2 = ^2(ei, e) ^i(ei)- In [|I^, the method has been apphed to deform the 
cychc supersymmetric chains, but we shall show that it can replace completely all 
other techniques in constructing quite arbitrary susy sequences. 

Indeed, the method can be now repeated by induction, leading to the sequence 
of first order intertwining operators 

Ak{e) = -^ + /3kix,e), fc = l,2,...,n, (14) 

where each (3^, k > 2, is simply the result of a finite difference operation performed 
on Pk-i- 

P,{x,e) = -A_i(a:,e,.i) - ^(g, 1 - e) 

f3k-i{x,ek-i) - /3fc-i(x,e) 

We have adopted here a shortcut notation making explicit only the dependence of 
Ak and f3k on the factorization constant introduced in the very last step, keeping im- 
plicit the dependence on the previous factorization constants (henceforth, the same 
criterion will be used for any other symbol depending on k factorization energies): 

Ak{e) = Ak{ei, ek-i, e), Pk{x, e) = pk{x, ei, e^-i, e). 

The Pk functions constructed at each next susy step automaticaly solve the Riccati 
equation with the potential Vk-i of the previous step 

- Pkix, e) + Plix, e) = 2[Vk.,{x, e,_i) - e], (16) 

or equivalently 

- piix, e) + Pl{x, e) = pi_^{x, ek-i) + Pl_i{x, ek-i) + 2{ek-i - e) (17) 
permitting to define the new potential 

Vfc(x, e) = Vfe„i(x, ek-i) + P'kix, e). (18) 
In the resulting sequence of the new Hamiltonians 

Hk{e) = '1^ + Vkix, e) = ^A,(e)4(e) + e, k = 1,2, n (19) 

each next is intertwined with the previous one, 

Hkie) Ak{e) = Afc(e) Hk^i{ek-i), Ho{eo) = Hq. (20) 

All of them are the higher order susy partners of the initial Hq 

Hkie)Ak = AkHo, (21) 

(compare with 0) where Ak is the higher order intertwining operator obtained in 
the factorized form: 

Ak = Ak{e)Ak-i{ek-i) ■ ■ ■ ^2(62)^1(61). (22) 
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(Here we follow our notation, making explicit the dependence of each Ak on its last 
factorization energy.) Once having the sequence of superpotentials /3fc(x, e), con- 



structed by the algorithm (p!5D , it is easy to determine the new eigenf unctions Uk{x, e) 
injected by the process into the Hamiltonians Hk-i{ek-i). Indeed, the Uk{x,ek) are 
defined by the first order differential equations 

Ak{e)ukix,e) = -^Ukix,e) + l3kix,e)ukix,e) = 0, /c = 1, 2, n. (23) 

Up to now we have maintained implicit the element of arbitrariness in the /3-functions 
at every recurrence step. However, the true advantages of the method can be seen 
if Pi{x,e) is wider than just a particular solution of (H). Indeed, for each fixed e, 
we have to our disposal the general solution of the Riccati equation (0), i.e., an 
entire one-parametric family of functions. Thus, e.g., the general solution of (§), for 
a given ei, depends on an additional integration parameter Ai. Given any particular 
solution f3ip[x,ei), the general one is obtained by an elementary transformation: 

/3i{x,ei) = /3ip{x,ei) - Xi - J exp ^^2 J f3ip{y , ti) dy^ dx = (3{x, ei, Xi). 

(24) 

Now, the algorithm (^) can be interpreted as an operation on two classes of func- 
tions Pi{x, e) and Pi{x, ei): an arbitrary member of the class Pi{x, e) subtracted from 
an arbitrary member of l3i{x, ei) in the denominator of ([TT|) . Quite similarly, the gen- 
eral iteration step (0) can be applied for /3fc_i(e) and Pk-i{^k-i) meaning any two su- 
perpotentials of the previous step, the only condition stating that they must solve the 
corresponding Riccati equation with the same potential Vk-2 (a fact which is auto- 
maticaly assured in the iterative process). If so, the method accumulates the integra- 
tion constants, permitting to recover the full parametric dependence on the n-th susy 
potentials in a purely algebraic way. Thus, e.g., [32{x, e) := [3{x;ei, e; Ai, A) with the 
dependence on ei, Ai, e, A inherited from the previous step (Ai and A selecting two so- 
lutions among the elements of the families /3i(x, ei) and (3i{x, e) respectively). In gen- 
eral, our algorithm yields l3k{x, e) := [3{x; ei, eu-i, e; Ai, Afc_i, A) where the pa- 
rameters e, A will be varied to perform the next step, whereas ei, eu-i, Ai, Xk-i 
will stay fixed. In what follows, we shall employ the symbols j3k{x, e) whenever want- 
ing to abbreviate, but we shall hold the dependence on the accumulated integration 
constants ei, e^-i, Ai, Afc_i if necessary. Let us now introduce the functions: 



Vtk{x,e) 



I3i{x, e), for k = 1; 

(25) 

^f3k{x,e) + (3k-i{x,ek-i), for A; > 2. 



Observe their fractal structure accumulating the finite difference derivatives (com- 
pare with [ p!3[] ) 

^k{x, e.) = , 2(e, ,-e,) ^ ^ 2. (26) 

\lk-i(x,ek-i) -\lk-i{x,ek) 
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Now, ( [T3| ) can be rewritten as Pk{x,e) = flk{x,e) — Pk-i{x,ek-i)- For example, for 
k = 2 we get (|lT]), while for /c = 3 we have 



l3{x; ei, €2, €3] Ai, A2, A3) =P{x; ei, Ai) + 



2(ei - es) 



f3{x; ei; Ai) - 62; X2) 
£2 - £3 



ei - 62 



ei - 63 



.(27) 



/?(x; ei; Ai) - (3{x; 62] X2) Pix; ei, Ai) - P{x; £3; A3) 



As one can also notice, the "n-susy potentials" Ki(x), present the different fractal 
structures for n even and odd. Summarizing: 



Vn{x, e„) = Vo{x) + J2 P'ki^, efc) = Vo{x) + E 

fc=l k=l 



COS 



(k + n)7r^ 



dx 



fifc(x,efc). (2J 



This expression has been already applied to simplify the study of ra-susy partners of 
the harmonic oscillator |[T^ and the Coulomb potential |19|. Notice now that when 
Afc = Afc_i, and taking the limit e^^i ^ e in (p!5D, the method produces a differential 
algorithm (a confluent susy operation) 



PrUx,e) = -Pk-iix,e) 



' d(3k-i{x,e) 
de 



(29) 



which can be freely iterated. Likewise, the finite differences can be replaced by the 
e-derivatives in one or more places of the fractal formulae (P7[). The mechanism 
has been implicitly applied by Stalhoffen (see Section III of [|^), in the context of 
Crum's determinant. 



3 Simple and confluent operations 

The difference-differential algorithms (^), (^) simplify remarkably the construction 
of transparent wells ^ In fact, the Riccati equation (|) for Vo(x) = has 
the general solution 

f3i{x;e;a) = -V2^ cot[V2^{x - a)], (30) 

with four different real 1-susy branches: S (singular), R (regular), P (periodic), and 
N (null), for different (real or complex) values of e and a (see Table I). For exam- 
ple, the regular case (R) leads to the modified Poschl- Teller potential V^(x, e, b) = 
—K^l cos\?[k(x + h)], while the null case (N) corresponds to the potential barrier 
Vf(x,e = 0,a) = (x — a)~^. All of them are particular cases of the Weierstrass 
p-function: 

Vi(x; e; a) = p(x + a, 3(4e/3)^, (4e/3)^) + ^. (31) 
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TABLE I. The four different real superpotentials /3i coming out from (^0|), 
depending on the values of e and the integration parameter a. In each case S 
means singular, R regular, P periodic, and N null. The parameters a and h are 
arbitrary real numbers. 



Case 


e 




a 


f3i{x,e) 


s 


e < 


e = in 


a 


—K coth[K(x — a)] 


R 


e < 


i^2\e\ = m 


b 


—K tanh[/t(x + b)] 


P 


e > 




a 


—k cot[A;(x — a)] 


N 








a 


1 



X — a 



By applying the algorithm (|TI]) to the superpotentials R and S as given in Table I, 
one immediately gets the transparent double wells in form of the quadratic Bargmann 
potentials 



V2(x,e2] 



(see Figure 1) 



[K^ K2) 



k1 csch^ [ki(x 



b)] + K2 sech^ [k2{x — a)] 



-Ki coth [ki{x + b)] + K2 tanh [k2{x — a)] 



(32) 




Figure 1: The double transparent wells (32) illustrate an application of the 
auto-Backlund transformation (pT|). The figure shows different members of the 
2-susy family ^2(2;, £2, ei) for ei — —4, 62 = —1, with integration parameters 
(a) a = 6 = 0, (b) a = 0.254, b = -1.018, and (c) a = 0.565, b = -2.262. 



The double wells (^2|) are well known in supersymmetry, but the use of the 
Backlund transformation makes the derivation considerably simpler, explaining also 
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their place in the theory of the KdV equation, vt — QvVx + Vxxx = 0. As a matter of 
fact, all the transparent n-susy wells can be now reinterpreted as the instantaneous 
forms of the multisoliton solutions propagating according to KdV (thus, e.g., the 
2-susy potentials ( p2l) originate the propagating soliton pairs [^). Presumably, 
the susy partners of an arbitrary Vo{x) generate as well the traveling "solitonic 
deformations" on a more general, non-null background. 

The operational advantages of the differential algorithm (p9D deserve as much 
attention. The application of the finite difference algorithm to the null wells is so 
easy, because for Vq = one knows the solutions of the Riccati eq. (^) for any e. For 
non-null Vq, ([TTl) can be applied efficiently only if one knows at least two solutions 
of (Q) for two different "factorization energies" e (compare [0, [|19|)- For a general 
Vq, however, this is not granted. It is therefore interesting, that the difficulty is 



circumvented by the differential algorithm (pQ]) which gives an easy access to higher 
order susy steps even if one starts from just one superpotential Pkix,e) for a fixed 
e. Indeed, by taking the limit e — > e^-i, one reduces (|1^ to the following equation 
for pk = Pkix,e = efc-i): 



-P'k + Pi = PLi + Pi 



l5 



k = 2,3, 



(33) 



which has the particular solution (3^ 

Pk = —Pk-1 + 



—Pk-1 and integrates easily yielding: 

g-2 / Pk-idx 



-A 



d 



dx 



In 



(34) 



Note that the right hand side of (|34|) depends essentially just on one integration 
constant. The recurrence (^) acquires even simpler form in terms of the new 'key 
functions' Bn = exp(— / Pndx + c„) (the integration constants Kn = expc„ are 
superfluous to determine /3„). In fact, integrating both sides of (0), taking to the 
exponent, differentiating again, and choosing properly the (inessential) constants 
Kn, one has: 

— (BnBn- 

dx 

When this method is applied to Vo{x) 
covers the Abraham-Moses oscillator |^ 
2-parametric family 



(35) 



: x^/2 with e = —1/2, the first step re- 
while the second step leads to the new 



Vo(x) 



X 

T 



d^ 

dx'^ 



In 



r, - 



dyey Ti- 



erf X 



(36) 



which are different from the potentials discussed in [|T^, because they appear for a 
fixed e and do not add 2 new spectral levels to Vq; we therefore propose to call them 
the 2-nd order Abraham-Moses potentials. The detailed analysis shows that (^) are 
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singular unless Ti = | and T2 > \ ^vr In 2, or Fi = — | y/n and r2 < —\ a/tt In 2. 
A nonsingular case with Ti = | a/tt and r2 = 0.308 is shown in Figure 2. The 
higher order Abraham- Moses functions can be as easily obtained in terms of (|3^) . 
Generically, our algorithm is close to the nice idea of Leble (though Leble did 
not construct the non-trivial contributions for 6 = e — e^^i = 0). 




Figure 2: The solid line is a nonsingular 2-susy confluent partner for the 
harmonic oscillator with factorization energy e — —1/2 (equation (|36|)), for 
Ti = V7r/2 and r2 = 0.308. The dashed line represents the harmonic oscillator 
potential 



An equally interesting application of the confluent formula ( PP]) arises for the 
periodic superpotential P in Table I. In this case, the first order susy partner is 
given by 

Vr\x, e) = V,ix, e) = (37) 

sm [k[x — a)\ 

while the corresponding second order susy step leads to 

V^^conf^^ ^) _ g^2 1 - COS [2k{x - g)] - kjx - g) sin [2k{x - a)] ^^^^ 



sin [2k{x — a)] — 2k{x — a) 



Notice that (|37| ) is a periodic potential with singularities at the points x„ = a+2n/k, 
n G Z, while V2°^^{x,e) has only one singularity. The confluent susy step (|29| ) has 
removed all the singularities from Vi°'^^{x,e) except that placed on x = g. The 
potential (|38|) has been plotted on Figure 3 (dashed curve) for g = 7, = 1/^2. 
Observe that we have rederived the Stahlhofen function, but the method is simpler 
(see equations (50) and (40) of ||2^). Moreover, some other potentials with analogous 



properties can be immediately obtained. Our Figure 3 reports the fourth order 
confluent potential V^°^^{x,e), with a = —7,k = l/\/2. Notice the appearance of 
two new symmetric singularities. 
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Figure 3: Examples of null, confluent n-susy potentials descending from the 
(P)-branch: V^°'^^{x) = V^°'^^{x;0.25;7), and V7°"f(x) = 1^4™"* (a;; 0.25; -7). 
For n even in (p8|), the number of singularities grows with n producing po- 
tentials of Stalhoffen type. The odd cases have always an infinite number of 
singularities. 



As becomes obvious, for the periodic superpotentials /3i in Table I, the method 
leads to two distinct classes of n-th order confluent potentials. For n even, the 
periodic superpotential Pi does not contribute directly to the sum (|28| ) (the term 
with k = 1 vanishes). The resulting susy partners have only a finite number of 
singularities. On the other hand, for n odd, the function f3i contributes to the first 
term in the sum ( PS]) and its global effect is never canceled. The corresponding susy 
partners Vi°'^^{x, e) therefore have infinite sequences of singularities. This difference 
reflects the distinct fractal structure of (^) for n even and odd. The possibility of 
mixed applications of ([15|) and (p9|) is open. 
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